On inverse problem of dynamics 

M. Rudnev*and V. Tent 

in 

§ . February 13, 2004 

(N 
X) 

' Abstract 

We study the question whether for a natural Hamiltonian system on a two-dimensional com- 
pact configuration manifold, a single trajectory of sufficiently high energy is almost surely enough 
, to reconstruct a real analytic potential. 

Q 

^ . AMS subject classification 37J05, 70H12 

Consider a compact configuration manifold M" (essentially we deal with the dimension n = 2) and 
a natural Hamiltonian system thereon, with Hamiltonian 

>■ H{p,q) = {p,p), + U{q), {p,q)eT*M''. 

(N 

■ Above, (■, ■)q is a Riemannian metric on M" and U a potential. The direct problem of dynamics on 

M"' is finding the trajectory q{t) C M", with initial conditions q{0) = qo and g(0) = vq, moving in 
I the known force field f{q) = —'VqU{q) on M", where the gradient Vg has been associated with the 

>n ■ metric (•, ■)q. 

Let us call the inverse problem of dynamics the problem of reconstruction of the potential by 
I observing the system's trajectories q{t). The first problem of this type was explored in Newton's 

a • Principia, in quest for a physical law determining the planetary motion compatible with observa- 

tional data^. In the general case, knowledge of infinitely many trajectories is required to completely 
solve the problem. In this note we show that in the special case when M" is two-dimensional, com- 
pact and topologically non-trivial, a single trajectory with sufficiently large energy would almost 
surely suffice to reconstruct the potential. 
I In the sequel, we assume that M" as well as all the quantities involved are real-analytic. Also 

suppose, there is an a-priori estimate < Co, ^q € M"', and we consider only the trajectories 

q{t) with total energy E > Cq, so M" coincides with the domain of possible motions. 

Theorem. Let n = 2, suppose is not diffeomorphic to S'^ or MP^. Almost every trajectory 
Q{t)i t ^ 0, with energy E > Cq, suffices to reconstruct the potential U as a real-analytic function 
on M^. 
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Let us recall the definition of a key set, or set of uniqueness, see e.g. 

Definition. Let D be a domain in M" and C^{D) the class of real-analytic functions in D. A 
set K d D is a key set if any f E C^{D) vanishing identically on K , vanishes identically on D. 

The theorem will follow from the following lemma. 

Lemma. If a real-analytic dynamical system A : x = F{x) on a compact phase space P, is 
non-singular (i.e. for no x P, F{x) = 0) and has a positive measure of closed orbits, then all its 
orbits are closed. 

Proof of the lemma 

Since the manifold P is compact and the vector field F is non-singular, there exists a Riemannian 
metric g on P, such that in this metric the vector field F has unit length at every x P. Fur- 
thermore, by compactness of P, the curvature (associated with g) of integral trajectories of A is 
bounded from above, and therefore there exists some > 0, such that any periodic orbit of A has 
period not smaller than T^. 

Let us partition the range [Tm,oo) of possible periods (henceforth periods stand for minimum 
periods) for closed orbits on intervals of some small length 6i to be specified. Let Ik = [Tm + 
k6i,Tm + {k + = [Tfe,Tfc_|_i), for k = 0, 1, ... . Let be the set of all closed orbits, whose 

periods lie in Ik- Then for some k = k^ the set = F^^, considered as a subset of P, has positive 
measure. (We refer to the sets F either as point sets or sets of orbits, depending on the context. As 
there are only measure-theoretical considerations involved, this should not cause confusion.) 

Let Ds^ix) be a codimension one disk in P, centered at some x G P, with radius 62 and 
perpendicular (in the sense of metric g) to the vector field F{x) at the point x. Let 82 be small 
enough, so that the vector field is transversal to (x) at every point of the disk. Clearly, 82 can 
be taken as a universal constant independent of x. 

Let xq be a Lebesgue point of the set F*. Recall that at a Lebesgue point, the density of the set 
is one. Let 70 be the closed trajectory passing through xq, so 70 € F^,. Take a disk D^^ixo). 

Then there is a well defined analytic Poincare map S from a disk D^{xq) contained in Ds^{xq), 
where e < ^Jj,^ 62, for some Ci = Ci{6i), and xq is a fixed point of S. On the disk D^{xo), points 
which are initial conditions for orbits from F^, form a set of positive measure, as xq is a Lebesgue 
point. Besides, the quantity 5i < Tm can be chosen small enough to ensure that any x G D£(xo)nF^, 
is also a fixed point of the map S. 

The union of all x G D^{xq) n F* is a positive measure subset of D^{xo), and hence is a key 
set (see [H] for the proof that every set of positive measure is a key set). Therefore, every point of 
De(xo) is an equilibrium of the map S, and hence an initial condition for a periodic orbit of A. Let 
Ff be the union of all such orbits, with initial conditions in D^(xo). Let F C P be the maximum 
connected open set, which contains F^ and is a union of periodic orbits. 

To complete the proof, let us show that the set F does not have a boundary, i.e. F = P. To 
show it, we use the following Gronwall type estimate. 

Proposition. Let (j){s) CT be a curve of length L, where s is a natural parameter with respect 
to the metric g. Then for all s G [0, L] and some absolute constant C2, 

T{<p{s)) < r(</.(0))e^2^, 

where T{<j)(s)) is the period of the closed orbit passing through the point (t){s). 
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Indeed, the proposition follows immediately from the following infinitesimal estimate: for some 



C2, 



< C2T{<t>{s)). 



Returning to the proof of the lemma, suppose the boundary dT is non-empty. As dV is a compact 
set, the distance (in the sense of metric g) between dV and the above mentioned point xq attains its 
minimum at some point y G dV. Connect xq and y by a geodesic segment. Let the latter segment 
have length L; clearly all its points, except y belong to T. Let 7^ be the trajectory of A with initial 
condition y. 

For any point xi 7^ y on the above geodesic segment, there is a uniform bound for the period of 
the corresponding closed orbit, by the proposition. Hence, for any such xi, there exists a uniform 
e (one can take e = ee~^^^ , where e has been defined earlier) such that an analytic Poincare map 
can be defined in exactly the same way as S above, but now with the domain De(xi). Choosing 
xi such that the intersection 7^ n Df,{x\) is not empty and repeating the key set argument leads 
to contradiction: all orbits in some tubular neighborhood of 7j^, including 7^^ itself, are periodic. 
Q.E.D. 



Proof of the theorem 

Consider a randomly chosen trajectory 7 on some energy level H~^{E), E > Cq, which is obviously 
a non-critical level. 

According to the lemma, either (i) all the trajectories on H~^{E) are periodic, or (ii) a randomly 
chosen initial condition {po,Qo) G H^^{E) results in a phase trajectory of infinite length almost 
surely. 

The former case (i) may occur only if is a so-called P-manifold. Indeed, according to the 
Maupertuis principle, the phase trajectories of motions with total energy E project onto as 
geodesic lines of the corresponding (Riemannian) Jacobi metric. P-manifolds are Riemannian man- 
ifolds, all whose geodesies are closed, see p. Topological properties of P-manifolds are characterized 
in great detail in various dimensions, within the framework of the Bott-Samelson theorem. In our 
(simplest possible) case, it is easy to see that can only be diffeomorphic to either 5^ or MP^. 

Indeed, the proof of the lemma implies that all the (closed) phase trajectories on the energy 
level E are homotopic to one another. Then their images on M"^ (under natural projection) are also 
homotopic. On the other hand, if is different from S'^ or RP^, the number of generators for its 
fundamental group equals at least two. As for any Riemannian metric there are closed geodesies in 
each free homotopy class, we would have a contradiction, unless is 5^ or MP^. 

In the case (ii), let q{t) = {qi{t) , q2{t)) be a randomly chosen trajectory: it almost surely has 
infinite length. Clearly, we can easily derive the gradient of U at every point of the trajectory from 
the Euler-Lagrange equations. 

Thus to complete the proof of the theorem, let us show that any non-closed trajectory is a 
key set. Consider orbit segments {q{t), t G [k,k + 1)}, k = 0, 1, . . . in M^. (Note that time can 
always be scaled to ensure that each segment is a simple curve in M^, or shorter time intervals can 
be considered.) As is compact, there is a limit point of the point sequence {q{k + 1/2)}. 
Consider a sufficiently small circle centered at q^:. There are two options. Either the circle intersects 
the trajectory q{t) at infinitely many points, or at some point on the circle, the trajectory q(t) 
intersects itself (transversely as it is a geodesic) infinitely many times. In the later case, take the 
point of infinite self-intersection for q^,, otherwise leave q^: as it is. In either case, there exists a point 
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g*, with the property that any sufficiently small circle centered at is intersected by the trajectory 
q{t) infinitely many times. 

Therefore, the force / and the potential U can be uniquely reconstructed on any sufficiently 
small circle centered at q.^ (an infinite point set on a circle is a key set for the circle), and therefore 
in some neighborhood of g*, and hence on the whole configuration manifold M^. Q.E.D. 

In conclusion, let us make several remarks. 

1. The theorem can restated in terms of analytic geodesic flows on compact Riemannian 2- 
manifolds. Namely, if is not diffeomorphic to a sphere or real projective plane, a randomly 
chosen geodesic suffices to reconstruct the metric, almost surely. Indeed, any Riemannian met- 
ric is locally conformal to the Euclidean one, i.e. can be locally associated with Hamiltonian 
H{p, q) = e^^''^''^^\p1 +P2)- Our theorem enables one to reconstruct the real-analytic function 
p{q1:Q2) locally near from the Hamilton equations, with subsequent analytic continuation 
to get the metric globally on M^. 

2. The theorem is essentially two-dimensional, as for n > 3 the fact that a random trajectory has 
infinite length does not suffice to reconstruct the potential. Consider for instance the Euler 
top, where = SO{3). In this case, the phase space is foliated by invariant two-tori, where 
the trajectories are in general conditionally periodic. Clearly, a projection of a single invariant 
two-torus onto the three-dimensional configuration space is not a key set. 

3. The lemma does not apply to the case of invariant tori of dimension higher than one. Indeed, 
a particular case of KAM theorem tells one that a sufficiently small perturbation of a non- 
degenerate Liouville-integrable Hamiltonian system in T*M^ yields a positive measure set of 
invariant n-tori, which do not fill the whole energy surface however. 

4. In the special case = RP^, the only possibility to have all the geodesic closed is the 
standard metric, by the theorem of L. Green ([!]). The case = 5^ has been a subject of 
extensive research for over 100 years, arguably beginning with the doctoral thesis of O. Zoll 
(|7j)- The reader is referred to the excellent book by L. Besse (jlj), which gives the issue a 
thorough treatment. 

5. The condition E > Cq also seems unavoidable, as for small energies the domain of possible 
motions can be a disk, with the Jacobi metric degenerate on the boundary, in which case one 
may expect the scenario, similar to the case = S^. 

6. Observe that in the exceptional case when is a P-manifold, the geodesic flow thereon is 
completely integrable (see e.g. for the proof of this fact). Hence our theorem implies that 
if n = 2, it is sufficient for restoration of the potential almost surely from a single trajectory 
that the system possess no other analytic integrals of motion but energy. It seems likely that 
in the latter weaker formulation the theorem should extend to the case n > 2, however we do 
not know how to prove it. 
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